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Abstract 

The reconstruction of the Higgs potential in the Minimal Supersymmetric Standard 
Model (MSSM) requires the measurement of the trilinear Higgs self-couplings. The 
'double Higgs production' subgroup has been investigating the possibility of detecting 
signatures of processes carrying a dependence on these vertices at the Large Hadron 
Collider (LHC) and future Linear Colliders (LCs). As reference reactions, we have 
chosen gg — > hh and e + e~ — > hhZ, respectively, where h is the lightest of the MSSM 
Higgs bosons. In both cases, the Hhh interaction is involved. For mu > 2mh, the 
two reactions are resonant in the H — > hh mode, providing cross sections which are 
detectable at both accelerators and strongly sensitive to the strength of the trilinear 
coupling involved. We explore this mass regime of the MSSM in the h — > bb decay 
channel, also accounting for irreducible background effects. 
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1 Introduction 



Considerable attention has been devoted to double Higgs boson production at future e + e~ 
and hadron colliders, both in the Standard Model (SM) and the MSSM (see Refs. [|I|, [| for 
an incomplete list of references for SM e + e~ and hadron colliders, respectively, and similarly 
H U for the MSSM). For the SM, detailed signal-to-background studies already exist for 
a LC environment ||, for both 'reducible' and 'irreducible' backgrounds [0, [7|, which have 
assessed the feasibility of experimental analyses. At the LHC, since here the typical SM 
signal cross sections are of the order of 10 fb [||, high integrated luminosities would be 
needed to generate a statistically large enough sample of double Higgs events. These would 
be further obscured by an overwhelming background, making their selection and analysis in 
a hadronic environment extremely difficult. Thus, in this contribution we will concentrate 
only on the case of the MSSM. 

In the Supersymmetric (SUSY) scenario, the phenomenological potential of these reac- 
tions is two-fold. Firstly, in some specific cases, they can furnish new discovery channels for 
Higgs bosons. Secondly, they are all dependent upon several triple Higgs self-couplings of 
the theory, which can then be tested by comparing theoretical predictions with experimental 
measurements. This is the first step in the reconstruction of the Higgs potential itselfQ. 

The Higgs Working Group (WG) has focused much of its attention in assessing the 
viability of these reactions at future TeV colliders. However, the number of such processes is 
very large both at the LHC and a LC [|[ ||, so only a few 'reference' reactions could be studied 
in the context of this Workshop. Work is in progress for the longer term, which aims to cover 
most of the double Higgs production and decay phenomenology at both accelerators 0. 

These reference reactions were chosen to be gg hh for the LHC (see top of Fig. [l]) 
and e + e~ — > hhZ for the LC (see bottom of Fig. p]), where h is the lightest of the MSSM 
scalar Higgs bosons. The reason for this preference is simple. Firstly, a stable upper limit 



exists on the value of m/,, of the order of 130 GeV, now at two-loop level |1(J, so that its 
detection is potentially well within the reach of both the LHC and a LC. In contrast, the 
mass of all other Higgs bosons of the MSSM may vary from the electroweak (EW) scale, 
0(m z ), up to the TeV region. In addition, as noted in Ref. ||, the multi-fe final state in 
gg — > hh — > bbbb, with two resonances and large transverse momenta, may be exploited in 
the search for the h scalar in the large tan /3 and moderate tua region. This is a corner of the 
MSSM parameter space that has so far eluded the scope of the standard Higgs production 



and decay modes JTTJ. (The symbol A here denotes the pseudoscalar Higgs boson of the 
MSSM, and we reserve the notation H for the heaviest scalar Higgs state of the model.) 
However, this paper will not investigate the LHC discovery potential in this mode, given the 
very sophisticated treatment of the background (well beyond the scope of this note) required 
by the assumption that no h scalar state has been previously discovered (see below). This 
will be done in Ref. ||. Furthermore, the gg — > hh and e + e~ — > hhZ modes largely 
dominate double Higgs production || |8|, at least for centre-of-mass (CM) energies of 14 TeV 
at the LHC and 500 GeV in the case of a LC, the default values of our analysis. (Notice 

lr The determination of the quartic self- interactions is also required, but appears out of reach for some 
time: see Refs. fa, 0] for some cross sections of triple Higgs production. 
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that we assume no polarization of the incoming beams in e + e~ scatterings.) Finally, when 
mjj > 2m/,, the two reactions are resonant, as they can both proceed via intermediate states 
involving H scalars, through gg —>■ H and e + e~ — > HZ, which in turn decay via H — > hh 
T^| . Thus, the production cross sections are largely enhanced H || (up to two orders of 



magnitude above typical SM rates at the LHC ||) and become clearly visible. This allows 
the possibility of probing the trilinear Hhh vertex at one or both these colliders. 

The dominant decay rate of the MSSM h scalar is into bb pairs, regardless of the value 



of tan (3 |12] . Therefore, the final signatures of our reference reactions always involve four 
6-quarks in the final state. (In the case of a LC environment, a further trigger on the 
accompanying Z boson can be exploited.) 

If one assumes very efficient tagging and high-purity sampling of 6-quarks, the background 
to hh events at the LHC is dominated by the irreducible QCD modes fll3| . Among these, 
we focus here on the cases qq, gg — > bbbb, as representative of ideal 6-tagging performances. 
These modes consist of a purely QCD contribution of 0{a^), an entirely EW process of 
C(«g m ) (with no double Higgs intermediate states) and an 0{a 2 s a 2 err ^ component consisting 
of EW and QCD interactions. (Note that in the EW case only qq initiated subprocesses are 
allowed at tree-level.) For a LC, the final state of the signal is bbbbZ, with the Z reconstructed 
from its decay products in some channel. Here, the EW background is of 0{a b em ) away from 
resonances (and, again, contains no more than one intermediate Higgs boson), whereas the 
EW/QCD background is proportional to (a 2 al m ). 

In general, EW backgrounds can be problematic due to the presence of Z vectors and 
single Higgs scalars yielding bb pairs, with the partons being typically at large transverse 
momenta and well separated. In contrast, the QCD backgrounds involve no heavy objects 
decaying to bb pairs and are dominated by the typical infrared (i.e., soft and collinear) QCD 
behaviour of the partons in the final state. However, they can yield large production rates 
because of the strong couplings. 

In this study, we investigate the interplay between the signal and background at both 
colliders, adopting detector as well as dedicated selection cuts. We carry out our analysis 
at both parton and hadron level. The plan of this note is as follows. The next Section 
details the procedure adopted in the numerical computation. Sect. ^ displays our results 
and contains our discussion. Finally, in the last section, we summarise our findings and 
consider possible future studies. 



2 Calculation 

For the parton level simulation, the double Higgs production process at the LHC, via gg fu- 
sion, has been simulated using the program of Ref. to generate the interaction gg — > hh, 
with the matrix elements (MEs) taken at leading-order (LO) for consistency with our treat- 
ment of the background. We then perform the two h — > bb decays to obtain the actual 
46-final state. For double Higgs production at a LC, we use a source code for the signal 
derived from that already used in Ref. M. At both colliders, amplitudes for background 



events were generated by means of MadGraph |T5j and the HELAS package [ffB . Note that 
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Figure 1: Diagrams contributing to gg — > hh (top) and e + e — > fo/iZ (bottom) in the 
MSSM. 



interferences between signal and backgrounds, and between the various background contri- 
butions themselves, have been neglected. This is a good approximation for the interferences 
involving the signal because of the very narrow width of the MSSM lightest Higgs boson. 
Similarly, the various background subprocesses have very different topologies, and one would 
expect their interferences to be small in general. 

The Higgs boson masses and couplings of the MSSM can be expressed at tree-level 
in terms of the mass of the pseudoscalar Higgs state, uia, and the ratio of the vacuum 
expectation values of the two neutral fields in the two iso-doublets, tan (3. At higher order 
however, top and stop loop-effects can become significant. Radiative corrections in the one- 
loop leading mf approximation are parameterised by 



3Gpmf 



log 



\/2tx 2 sin 2 (3 rn 



(1) 



where the SUSY breaking scale is given by the common squark mass, m s , set equal to 1 TeV 
in the numerical analysis. If stop mixing effects are modest at the SUSY scale, they can 
be accounted for by shifting m| in e by the amount Am| = A 2 [l — A 2 /(12m 2 s )] (A is the 
trilinear common coupling). The charged and neutral CP-even Higgs boson masses, and the 
Higgs mixing angle a are given in this approximation by the relations: 
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as a function of rriA and tan (3. The triple Higgs self-couplings of the MSSM can be param- 



eterised |L7], [18| in units of M z /v, 



246 GeV, as, 
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Next-to- leading order (NLO) effects are certainly dominant, though the next-to-next- 
to-leading order (NNLO) ones cannot entirely be neglected (especially in the Higgs mass 
relations). Thus, in the numerical analysis, the complete one-loop and the leading two- loop 
corrections to the MSSM Higgs masses and the triple Higgs self- couplings are included. The 
Higgs masses, widths and self-couplings have been computed using the HDECAY program 
described in Ref. which uses a running 6-mass in evaluating the h — > bb decay fraction. 
Thus, for consistency, we have evolved the value of m& entering the hbb Yukawa couplings of 
the h ^ bb decay currents of our processes in the same way. 

For our analysis, we have considered tan/3 = 3 and 50. For the LHC, high values of tan/? 
produce a signal cross section much larger than the tan (3 = 3 scenario, over almost the entire 
range of m^. However, this enhancement is due to the increase of the down-type quark-Higgs 
coupling, which is proportional to tan (3 itself, and serves only to magnify the dominance of 
the quark box diagrams of Fig. |I]. Unfortunately, these graphs have no dependence on either 
of the two triple Higgs self-couplings entering the gluon-gluon process considered here, i.e., 
Xhhh and \nhh- Thus, although the cross section is comfortably observable, all sensitivity to 
such vertices is lost. Therefore, the measurement of the triple Higgs self-coupling, Xnhh, is 
only feasible at the LHC for low tan/3 due to the resonant production of the heavy Higgs 
boson (see Fig. 5a of Ref. 0). 

In contrast, the cross section for double Higgs production at the LC is small for large 
tan (3 because there is no heavy Higgs resonance (see Fig. 8 of Ref. As soon as it becomes 
kinematically possible to decay the heavy Higgs into a light Higgs pair, the ZZH coupling 
is already too small to generate a sizable cross section. Furthermore, the continuum MSSM 
cross section is suppressed with respect to the SM cross section since the MSSM couplings 
ZZH and ZZh vary with cos(/3 — a) and sin(/3 — a), respectively, with respect to the 
corresponding SM coupling. Notice that in this regime, at a LC, the Xhhh vertex could in 
principle be accessible instead, since Xhhh 3> Xnhh (see Fig. 2 of Ref. 0) and because of the 
kinematic enhancement induced by rrih Unfortunately, we will see that the size of 

the e + e~ — > hhZ cross section itself is prohibitively small. 
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We assume that 6-jets are distinguishable from light-quark and gluon jets and no efficiency 
to tag the four 6-quarks is included in our parton level results. We further neglect considering 
the possibility of the 6-jet charge determination. Also, to simplify the calculations, the Z 
boson appearing in the final state of the LC process is treated as on-shell and no branching 
ratio (BR) is applied to quantify its possible decays. In practice, one may assume that 
it decays leptonically (i.e., Z — > with £ = e,fi,r) or hadronically into light-quark 

jets (i.e., Z — > qq, with q ^ b), in order to avoid problems with 66-quark combinatorics. 
Furthermore, in the LC analysis, we have not simulated the effects of Initial State Radiation 
(ISR), beamstrahlung or Linac energy spread. Indeed, we expect them to affect signal and 
backgrounds rather similarly, so we can neglect them for the time being. Indeed, since a 
detailed phenomenological study, including both hadronisation and detector effects, already 
exists for the case of double Higgs-strahlung in e + e~ 0, whose conclusions basically support 
those attained in the theoretical study of Ref . , we limit ourselves here to update the latter 
to the case of the MSSM. 



So far only resonant production gg — > H — > hh — > bbbb has been investigated [13], with 



full hadronic and detector simulation and considering also the (large) QCD backgrounds, 
and a similar study does not exist for continuum double Higgs production at the LHC. (See 
Ref. [[H| for a detailed account of the gg — > H — > hh —>■ 7766 decay channel.) The event 
simulation has been performed by using a special version of PYTHIA f2l| , in which the 



relevant LO MEs for double Higgs production of Ref. [14j have been implemented by M. El 
Kacimi and R. Lafaye. These MEs take into account both continuum and resonant double 
Higgs boson production and their interferences. (The insertion of those for e + e~ processes is 
in progress.) The PYTHIA interface to HDECAY has been exploited in order to generate 
the MSSM Higgs mass spectrum and the relevant Higgs BRs, thus maintaining consistency 
with the parton level approach. As for the LHC detector simulation, the fast simulation 
package was used, with high luminosity (i.e., J Cdt = 100 fb _1 ) parameters. 

The motivation for our study is twofold. On the one hand, to complement the studies 
of Ref. Jl3[ by also considering the continuum production gg — > hh — > bbbb at large tan /3. 



On the other hand, to explore the possibility of further kinematic suppression of the various 
irreducible backgrounds to the resonant channel at small tan /3. 



3 Results 

3.1 The LHC analysis 

In our LHC analysis, following the discussion in Sect. |2], we focus most of our attention on 
the case tan/3 = 3, with ttia = 210 GeV, although other combinations of these two MSSM 
parameters will also be considered. We further set A = —\i = 1 TeV and take all sparticle 
masses (and other SUSY scales) to be as large as 1 TeV. 
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3.1.1 gg — > hh — > bbbb at parton level 

In our parton level analysis, we identify jets with the partons from which they originate 
(without smearing the momenta) and apply all cuts directly to the partons. We mimic the 
finite coverage of the LHC detectors by imposing a transverse momentum threshold on each 
of the four 6-jets, 

p T (b) > 30 GeV (4) 
and requiring their pseudorapidity to be 

|tK&)|<2.5. (5) 

Also, to allow for their detection as separate objects, we impose an isolation criterium among 
6-jets, 

AR(bb) > 0.4, (6) 

by means of the usual cone variable AR(ij) = a/ ' Arj{ij) 2 + A(f>(ij) 2 , defined in terms of 
relative differences in pseudo-rapidity rjij and azimuth 0^ of the i-th and j-th 6-jets. 

As preliminary and very basic selection cuts (also to help the stability of the numerical 
integration), we have required that the invariant mass of the entire 46-system is at least 
twice the mass of the lightest MSSM Higgs boson (apart from mass resolution and gluon 
emission effects), e.g., 

mibbbb) > 2m h - 40 GeV, (7) 
and that exactly two /i-resonances are reconstructed, such that 

\m(bb) -m h \ < 20 GeV. (8) 

In doing so, we implicitly assume that the h scalar boson has already been discovered and 
its mass measured through some other channel, as we have already intimated in the Intro- 
duction. 

After the above cuts have been implemented, we have found that the two 46-backgrounds 
proceeding through EW interactions are negligible compared to the pure QCD process. In 
fact, the constraints described in eqs. produce the strongest suppression, almost 

completely washing out the relatively enhancing effects that the cuts in (@)-(@) have on 
the EW components of the backgrounds with respect to the pure QCD one, owning to the 
intermediate production of massive Z bosons in the former. In the end, the production 
rates of the three subprocesses scale approximately as their coupling strengths: i.e., 0{a^) : 
0{a 2 s a 2 ern ) : 0(a^ m ). Therefore, in the reminder of our analysis, we will neglect EW effects, 
as they represent not more than a 10% correction to the QCD rates, which are in turn 
affected by much larger QCD f^-factors. As for the pure QCD background itself, it hugely 
overwhelms the double Higgs signal at this stage. The cross section of the former is about 
7.85 pb, whereas that of the latter is approximately 0.16 pb. 

To appreciate the dominance of the rrih cuts, one may refer to Fig. ^ where the distri- 
butions in transverse momentum of the four p^-ordered 6-quarks (such that pr(^i) > ••• > 
Pribi)) of both signal and QCD background are shown. Having asked the four 6-jets of the 
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Figure 2: Distributions in transverse momentum of the four p^-ordered fe-jets in gg — > hh —> 
bbbb and in the QCD background, after the cuts (|)-@ at the LHC, for tan/5 = 3, m h = 104 
GeV and m# ~ 220 GeV. Normalisation is to unity. 

background to closely emulate the gg —>■ hh — > bbbb kinematics, it is not surprising to see 
a 'degeneracy' in the shape of all spectra. Clearly, no further background suppression can 
be gained by increasing the prip) cuts. The same can be said for 77(6) and AR(bb). Others 
quantities ought to be exploited. 

In Fig. ^, we present the signal and QCD background distributions in the minimum angle 
formed between the two 6-quarks coming from the 'same Higgs' (i.e., those fulfilling the cuts 
in (|5D) in the 46-system rest frame (the plot is rater similar for the maximum angle, thus 
also on average). There, one can see a strong tendency of the two 26-pairs produced in the 
Higgs decays to lie back-to-back, reflecting the 2 — > 2 intermediate dynamics of Higgs pair 
production via gg — > hh. Missing such kinematically constrained virtual state, the QCD 
background shows a much larger angular spread towards small 6 m i n {bb) values, eventually 
tamed by the isolation cut @. 

The somewhat peculiar shape of the signal distribution is due to destructive interference. 
Recall that the signal contains not only diagrams proceeding via a heavy Higgs resonance (the 
upper-left hand graph of Fig. [I]), which results in the large peak in Fig. |3|, but also contains 
a continuum contribution mediated by box graphs (the upper-right hand graph of Fig. [I]). 
These two contributions destructively interfere leading to the depletion of events between 
the large back-to-back peak and the small remaining 'bump' of the continuum contribution 
as seen in Fig. || 
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Figure 3: Distributions in minimum relative angle (in radians) in the 46-system rest frame 
between two 6-jets reconstructing in gg — > hh — > bbbb and in the QCD background, 
after the cuts ®-(D at the LHC, for tan (3 = 3, m h = 104 GeV and m H ~ 220 GeV. 
Normalisation is to unity. 

In the end, a good criterium to enhance the signal-to-background ratio (S/B) is to require, 
e.g., 9{bb) > 2.4 radians, i.e., a separation between the 26-jets reconstructing the lightest 
Higgs boson mass of about 140 degrees in angle. (Incidentally, we also have investigated the 
angle that each of these 26-pairs form with the beam axis, but found no significant difference 
between signal and QCD background). 

An additional consequence that one should expect from the presence of two intermediate 
massive objects in gg —>■ hh — > bbbb events is the spherical appearance of the jets in the 
final state, in contrast to the usual planar behaviour of the infrared QCD interactions. 
These phenomena can be appreciated in Fig. Notice there the strong tendency of the 
background to yield high thrust configurations, again controlled by the separation cuts when 
T approaches unity. On the contrary, the average value of the thrust in the signal is much 
lower, being the effect of accidental pairings of 'wrong' 26-pairs (the shoulder at high thrust 
values) marginal. An effective selection cut seems to be, e.g., T < 0.85. 

Furthermore, if the heavy Higgs mass is sufficiently well measured at the LHC then one 
can exploit the large fraction || of 46-events which peak at itlh in the signal, as dictated 
by the H — > hh decay, improving the signal-to-background ratio. This peak at m# can be 
clearly seen in the left hand plot of Fig. || where it dominates the QCD background, even for 
bins 13 GeV wide. In fact, not only could the QCD background be considerably suppressed 
but also those contributions to gg — > hh not proceeding through an intermediate H state 
should be removed, this greatly enhancing the sensitivity of the signal process to the Xnhh 
coupling. This can be seen in the right hand plot of Fig. [5] where the signal is shown on a 
logarithmic scale. The continuum contribution due to the box graphs (and its destructive 
interference with the heavy Higgs decay contribution) is now evident although one should 
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Figure 4: Distributions in thrust in the rest frame of the 46-system in gg — > hh — > bbbb and 
in the QCD background, after the cuts (H)-© at the LHC, for tan (3 = 3, m h = 104 GeV 
and tuh — 220 GeV. Normalisation is to unity. 
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Figure 5: Distributions in invariant mass of the 4fe-system in gg — > hh — > bbbb and in the 
QCD background, after the cuts @-@ at the LHC, for tan/3 = 3, = 104 GeV and 
mil — 220 GeV. Normalisation is to unity. The left hand plot shows both the signal (solid 
curve) and the QCD background (dashed curve), distributed in 5 GeV bins. The same signal 
is also shown as a histogram for a more experimentally realistic binning of 13 GeV. The right 
hand plot also shows the signal (collected in 5 GeV wide bins) on a logarithmic scale. Here 
the structure of the continuum contribution (and its destructive interference with the heavy 
Higgs decay contribution) can be seen. 

note that it is considerably suppressed compared to the peak at m#. 

Now, if a less than 10% mass resolution can be achieved on the light and heavy Higgs 
masses, then one can tighten cut @ to \m{bb) —m h \ < 10 GeV and introduce the additional 
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cut \m(bbbb) — m#| < 20 GeV. These cuts taken together with those in 9(bb) and T already 
suggested, reduce the QCD background to the same level as the signal. In fact, we have 
found that the cross section of the background drops to approximately 174 fb whereas that 
of the signal remains as large as 126 fb, this yielding a very high statistical significance at 
high luminosity Even for less optimistic mass resolutions the signal-to-background ratio 
is still significantly large. For example, selecting events with \m(bb) — rrih\ < 20 GeV and 
\m(bbbb) — w,h\ < 40 GeV, the corresponding numbers are approximately 102 fb for the 
signal and 453 fb for the background. Notice that the signal actually decreases as these 
Higgs mass windows are made larger. This is due to our insistence that exactly two bb pairs 
should reconstruct the light Higgs mass. As the light Higgs mass window is enlarged from 
rrih ± 10 GeV to ± 20 GeV, it becomes more likely that accidental pairings reconstruct 
the light Higgs boson. Since one is then unable to unambiguously assign the b quarks to the 
light Higgs bosons, the event is rejected and the signal drops. 

Although we have discussed here an ideal situation which is difficult to match with more 
sophisticated hadronic and detector simulations, it still demonstrates that the measurement 
of the Xuhh coupling could be well within the potential of the LHC, at least for our particular 
choice of MSSM parameters. Comforted by such a conclusion, we now move on to more 
realistic studies. 



3.1.2 gg — > hh — > bbbb at the LHC experiments 

Although the LHC experiments will be the first where one can attempt to measure the Higgs 
self-couplings, the analysis is very challenging because of the smallness of the production cross 
sections. Even in the most favourable cases, the production rate is never larger than a few 



picobarns, already including one-loop QCD corrections, as computed in Ref. [pf] . The cross 
sections at this accuracy are given in Tab. [T], for the resonant process (case 1 with mu = 220 
GeV) as well as three non resonant scenarios: one at the same tan f3 but with the H —>■ hh 
decay channel closed (case 2), a second at very large tan/? and no visible resonance (case 3) 
and, finally, the SM option (case 4, where mu identifies with the mass of the standard Higgs 
state). 



case 


model 


tan j3 


m h (GeV) 


A (TeV) 


H (TeV) 


a (fb) 


dominant mode 


1 


MSSM 


3 


104 


+ 1 


-1 


2000 


gg — > H —>■ hh 


2 


MSSM 


3 


100 


+1 


-1 


20 


gg -> hh 


3 


MSSM 


50 


105 


+1 


+1 


5000 


gg -> hh 


4 


SM 




105 






40 


gg -> hh 



Table 1: Cross sections for double Higgs production hh at the LHC via gluon-gluon fusion 
at NLO accuracy, for three possible configurations of the MSSM and in the SM as well. 
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3.1.3 LHC trigger acceptance 

For 46-final states, possible LHC triggers are high px electron/muons and jets. As an exam- 
ple, the foreseen ATLAS level 1 trigger thresholds on px and acceptance for a 46-selection 
(with the four 6-jets reconstructed in the detector) are given in Tab. |2], assuming the LHC 
to be running at high luminosity. 



trigger type: 


1 e 


1/i 


2 // 


1 jet 


3 jets 


4 jets 


total 


Pt in GeV 


30 


20 


10 


290 


130 


90 




case 1, e(bbbb) in % 


0.01 


0.01 


0.4 


0.08 


0.08 


0.05 


0.53 


case 2 


< 0.01 


< 0.01 


2.1 


2.9 


3.8 


4.2 


8.8 


case 3 


< 0.01 


< 0.01 


2.2 


2.7 


3.8 


4.1 


8.7 


case 4 


< 0.01 


< 0.01 


2.0 


2.5 


3.3 


3.6 


7.8 



Table 2: Kinematical acceptance of the ATLAS detector to trigger four 6-jets (including 
detector acceptance) at high luminosity. 

The overall trigger acceptance is at best 8-9%, for cases 2,3,4. The very low efficiency for 
case 1 is clearly a consequence of the small value of the difference tuh — ^rrih, translating into 
a softer prib) spectrum with respect to the other cases (compare the left-hand with the right- 
hand side of Fig. ^j. One can further see in the left-hand plot of Fig. [| that the bulk of the 
signal lies below the lowest pr(b) threshold of Tab. (i.e., 90 GeV), so that adopting smaller 
trigger thresholds could result in a dramatic enhancement of our efficiency. Of course, this 
would also substantially increase the low transverse momentum QCD background, as we can 
see in the parton level analysis of Fig. |2|. 

For example, by lowering the thresholds to 180, 80 and 50 GeV for 1, 3 and 4 jets, 
respectively (compare to Tab. ||), the overall trigger acceptance on the signal goes up to 
1.8%, i.e., by almost a factor of 4. Meanwhile though, the ATLAS level-1 jet trigger rates 
increase by a factor of 10 p2| . Anyhow, even for our poor default value of e(bbbb) in Tab. ^ 
we will see that case 1 still yields a reasonable number of events in the end. Optimisations 
of the 6-jet transverse momentum thresholds are in progress [Q. 

3.1.4 LHC events selection for gg — > hh — > bbbb 

Jets are reconstructed merging tracks inside AR(bb) = 0.4. Only jets with transverse en- 
ergy/momentum greater than 30 GeV and with \rj(b)\ < 2.5 are kept. (Thus, the same 
cuts as in the parton level analysis, now applied instead to jets.) The effect from pile up is 
included in the resolution. A jet energy correction is then applied. 

The invariant masses of each jet pair can then be computed. Assuming that the lightest 
Higgs boson mass is known, events with m(bb) sufficiently close to can efficiently be 
selected, see Fig. |7|. Another cut on the AR(bb, bb) between pairs of 6-jets can also be 
applied to reduce the intrinsic combinatorial background, since the latter concentrates at 
large AR(bb, bb) values, see Fig. ||. 
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Figure 6: Reconstructed transverse energy /momentum for 6-jets in gg — > hh — > bbbb events 
of case 1 (left plot) and 6-jets in gg — > hh — > bbbb events of case 2 (right plot) with ATLAS 
fast simulation [T2^] at high luminosity. Normalisation is arbitrary. 
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Figure 7: Reconstructed invariant mass distribution of 26-jet pairs in continuum gg — > 
hh — > bbbb events (case 2) with the fast simulation at high luminosity. Normalisation is 
arbitrary. (Results of a Gaussian fit are also given.) 

For already discussed, we can further impose that the invariant mass of the 

four 6-jets should be the heavy Higgs mass, m#, in order to select the H — > hh resonance, 
as confirmed by Fig. ^. In the other three cases, where the H —* hh splitting is no longer 
dominant (MSSM) or non-existent (SM), one can still insist that the 46-jet invariant mass 
should be higher than two times the lightest Higgs mass, see Fig. [H] and recall eq. (0). 
Finally, following Fig. [11], by constraining the 6-jets pairs four-momenta around the known 
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Figure 8: Reconstructed AR(bb, bb) between 26-jet systems from h — > bb decays in continuum 
gg — > hh — > bbbb events (case 2) with the fast simulation at high luminosity The dashed 
histogram shows the same distribution for all pairs of jets. Normalisation is arbitrary. 

light Higgs mass value, m/,, one can further reject the intrinsic background by means of the 
m(bbbb) spectrum. 

3.1.5 LHC 6-tagging in gg — > hh — > bbbb 

The 6-tagging efficiency at high luminosity is set to 50%, with dependent correction 
factors for jets rejection. An average rejection of 10 for c-jets and 100 for light-jets can be 
expected. We then studied the effect on the selection efficiency of requiring from one to four 
6-tags, although it is clear that, according to the parton level studies, the huge background 
rate demands four 6-tags, leading to a 6% tagging efficiency overall. 

3.1.6 Event rates at the LHC 

Taking into account all the efficiencies described above, and using the NLO normalisation of 
Tab. U, one can extract the number of expected events per year at the LHC at high luminosity 
given in Tab. The selection cuts enforced here are the following. For a start, we have kept 
configurations where \m(bb) — rrih\ < 30 GeV (cases 1,3,4) or \m{bb) — rrih\ < 20 GeV (case 2) 
and AR(bb, bb) < 2.5 (all four cases). (If more than two m^'s are reconstructed, the best two 
26-pairs are selected according to the minimum value of 5M 2 = [m^— m(bb)] 2 + [rrih— m'(bb) 2 ].) 
Then, a cut on m(bbbb) is applied: in presence of the H —> hh resonance (case 1) we have 
kept events within an m H mass window of ±2a (about 82% of the total number survive); 
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Figure 9: Reconstructed 46-jet invariant mass for 6-jets coming from the hh pair in gg — > 
hh — > bbbb events (case 1) with the fast simulation at high luminosity. The dashed histogram 
shows the same distribution for all groups of four jets. Normalisation is arbitrary. (Results 
of a Gaussian fit to the first spectrum are also given.) 

otherwise (cases 2,3,4) we have adjusted the m{bbbb) > 2m h cut so to keep 90% of the sample. 
In the end, one finds the numbers in Tab. ^ that are encouraging indeed. 





case 1 


2 


3 


4 


a in fb 


2000 


20 


5000 


40 


trigger threshold acceptance 


0.53% 


8.8% 


8.7% 


7.8% 


mass windows 


60% 


50% 


40% 


40% 


46-tagging 


6% 


6% 


6% 


6% 


events/year (no tagging) 


636 


88 


17400 


125 


events/year (four 6-tags) 


38 


5.3 


1044 


7.5 



Table 3: Total rates for gg — > hh — > bbbb, after all efficiencies have been included and 
selection cuts (§)-© enforced at hadron level, with 100 fb -1 per year of luminosity. 

In conclusion then, looking at the results in Tab. and bearing in mind the potential seen 
in reducing the pure QCD background via gg — > G{afy — > bbbb (see Figs. §]-§), one should be 
confident in the LHC having the potential to measure the \Hhh coupling in resonant H — > hh 
events (case 1). To give more substance to such a claim, we have now initiated background 
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Figure 10: Reconstructed 46-jet invariant mass for 6-jets coming from the hh pair in gg — > 
hh — > bbbb events (case 4) with the fast simulation at high luminosity. The dashed histogram 
shows the same distribution for all groups of four jets. Normalisation is arbitrary. 

studies at hadron and detector level, following the guidelines obtained by the parton level 
analysis |§. As for other configurations of the MSSM (such as case 2) or in the SM (case 4), 
the expectations are more pessimistic. Case 3 deserves further attention. In fact, notice the 
large number of events surviving and recall what mentioned in the Introduction concerning 
the potential of the non-resonant gg — *> hh — > bbbb process as a discovery channel of the light 
Higgs boson of the MSSM in the large tan (3 region at moderate values, a corner of the 
parameter space where the h coverage is given only by SM-like production/decay modes, 
thus not allowing one to access information on the MSSM parameters. Results on this topic 
too will be presented in Ref. |J. 

3.2 The LC analysis 

Here, we closely follow the selection procedure advocated in Ref. ]7|. In order to resolve the 
four 6-jets as four separate systems inside the LC detector region, we impose the following 
cuts. First, that the energy of each 6-jet is above a minimum threshold, 

E(b) > 10 GeV. (9) 

Second, that any 6-jet is isolated from all others, by requiring a minimum angular separation, 

cos 0(6, b) < 0.95. (10) 




16 



1600 



1400 



1200 



1000 



800 



600 



400 



200 




100 200 300 400 500 600 700 800 900 1000 
m(bbbb) distribution m ( GeV ) 



Figure 11: Reconstructed 46-jet invariant mass for 6-jets coming from the hh pair in gg — > 
hh — > bbbb events (case 4) with the fast simulation at high luminosity. Here, the energy of 
the jet pairs is recalculated using the m/, constraint. The dashed histogram shows the same 
distribution for all groups of four jets. Normalisation is arbitrary. 

Similarly to the hadronic analysis, one can optimise S/B by imposing the constraints 0, 

m{bbbb) > 2m h - 10 GeV, (11) 

\m(bb) -m h \ < 5 GeV, (12) 

on exactly two combinations of 26-jets. Here, note that the mass resolution adopted for the 
quark systems is significantly better than in the LHC case, due to the cleanliness of the e + e~ 
environment and the expected performance of the LC detectors in jet momentum and angle 



reconstruction EM. Thus, given such high mass resolution power from the LC detection 
apparatus, one may further discriminate between h and Z mass peaks by requiring that 
none of the 26-jet pairs falls around rriz, 



\m{bb) — mz\ > 5 GeV. 



(13) 



Moreover, in the double Higgs-strahlung process e + e _ — > hhZ, the four 6-quarks are pro- 
duced centrally, whereas this is generally not the case for the background (see the discussion 
in Ref. 0). This can be exploited by enforcing 



cos 6(bb, bbb, bbbb) | < 0.75, 



(14) 
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Figure 12: Cross sections in femtobarns for the e + e~ — > hhZ signal in the h — > bbbb decay 
channel, at a LC with 500 GeV as CM energy, as a function of for tan/3 = 3 and 
50. Our acceptance cuts in energy and separation of the four 6-quarks (P)-(|TU) have been 
implemented. No beam polarisation is included. 



where 6(bb, bbb, bbbb) are the polar angles of all two-, three- and four-jet systems. 

Fig. IT2] shows the production and decay rates of the signal process, e + e~ — > hhZ — > bbbbZ, 



as obtained at the partonic level, after the cuts (|9|)-(|T0D have been implemented. The MSSM 
setup here includes some mixing, having adopted A = 2.4 TeV and \i — 1 TeV, at both 
tan/3 = 3 and 50. Notice the onset of the H — > hh — > bbbb decay sequence in the Higgs- 
strahlung process e + e~ — > HZ at low tan/3. The same does not occur for large values, as 
previously explained. The impact of the above jet selection cuts on the signal is marginal, 
as the 6-quarks are here naturally isolated and energetic, being the decay products of heavy 
objects. In fact, the rates in Fig. |T2] would only be 10-20% higher if all the 46-quark phase 
space was allowed (the suppression being larger for smaller Higgs masses). At the height 
of the resonant peak around rrih ~ 104 GeV at tan (3 = 3, the signal rate is not large but 
observable, yielding more than one event every 1 fb _1 of data. For a high luminosity 500 
GeV TESLA design ||25|| , this would correspond to more than 300 events per year. Given 



the very high efficiency expected in tagging 6-quark jets, estimated at 90% for each pairs of 
heavy quarks p6| , one should expect a strong sensitivity to the triple Higgs self-coupling. 
The situation at large tan (3 is much more difficult instead, being the production rates smaller 
by about a factor of 10. 

In the left-hand side of Fig. ITS] we present the EW background, after the constraints in 



(p|)— fllQ|) have been enforced, in the form of the four dominant EW sub-processes. These four 
channels are the following. 

1. e + e~ — > ZZZ — > bbbbZ, first from the left in the second row of topologies in Fig. 3 of 
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Figure 13: Cross sections in femtobarns for the dominant components of the EW (left) and 
EW/QCD (right) background to the e + e~ — > hhZ signal in the h — > bbbb decay channel, 
at a LC with 500 GeV as CM energy, as a function of nth for tan (3 = 3 (top) and 50 
(bottom). Our acceptance cuts in energy and separation of the four 6-quarks (P) (pTOj) have 
been implemented. No beam polarisation is included. 

Ref. 0. That is, triple Z production with no Higgs boson involved. 

2. e + e~ — > h/HZZ — > bbbbZ, first(first) from the left(right) in the fifth(fourth) row of 
topologies in Fig. 2 of Ref. |7|] (also including the diagrams in which the on-shell Z is 
connected to the electron-positron line). That is, single Higgs-strahlung production in 
association with an additional Z, with the Higgs decaying to bb. The cross sections of 
these two channels are obviously identical. 

3. e + e~ — > h/HZ — > Z*Z*Z — > bbbbZ, first from the right in the third row of topologies 
in Fig. 2 of Ref. @ . That is, single Higgs-strahlung production with the Higgs decaying 
to bbbb via two off-shell Z* bosons. 

4. e + e~ — > Zh/H — ► bbZ*Z — > bbbbZ, first(first) from the right(left) in the first(second) 
row of topologies in Fig. 2 of Ref. . That is, two single Higgs-strahlung production 
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channels with the Higgs decaying to bbZ via one off-shell Z* boson. Also the cross 
sections of these two channels are identical to each other, as in 2. 

The 0{a 2 s al m ) EW/QCD background is dominated by e + e~ — ► ZZ production with one 
of the two Z bosons decaying hadronically into four fe-jets. This subprocess corresponds to 
the topology in the middle of the first row of diagrams in Fig. 4 of Ref. 0. Notice that 
Higgs graphs are involved in this process as well (bottom-right topology in the mentioned 
figure of [0). These correspond to single Higgs-strahlung production with the Higgs scalar 
subsequently decaying into bbbb via an off-shell gluon. Their contribution is not entirely 
negligible, owing to the large ZH production rates, as can be seen in the right-hand side of 
Fig. |13[ The interferences among non-Higgs and Higgs terms are always negligible. 

In performing the signal-to-background analysis, we have chosen two representative points 
only, identified by the two following combinations: (i) tan (3 = 3 and rriA = 210 GeV 
(yielding m h « 104 GeV and m H « 220 GeV); (ii) tan/5 = 50 and m A = 130 GeV (yielding 
rrih ~ 120 GeV and mn ~ 130 GeV). These correspond to the two asterisks in Fig. [L2], 
that is, the maxima of the signal cross sections at both tan f3 values. The first corresponds 
to resonant H — > hh production, whereas the latter to the continuum case. If we enforce 
the constraints of eq. (|TTD — , the suppression of both EW and EW/QCD is enormous, so 
that the corresponding cross sections are of O(10~ 3 ) fb, while the signal rates only decrease 
by a factor of four at most. This is the same situation that was seen for the SM case in 
Ref. f7|. Indeed, in the end it is just a matter of how many signal events survive, the sum 
of the backgrounds representing no more than a 10% correction (see Fig. 11 of Ref. 0). 
For example, after 500 fb -1 of data collected, one is left with 156 and 15 events for case (i) 
and (ii), respectively. However, these numbers do not yet include 6-tagging efficiency and Z 
decay rates. 

4 Summary 

To summarise, the 'double Higgs production' subgroup has contributed to the activity of 
the Higgs WG by assessing the feasibility of measurements of triple Higgs self-couplings at 
future TeV colliders. The machines considered were the LHC at CERN (14 TeV) and a 
future LC running; at 500 GeV. In both high luminosity setup was assumed, given 

the smallness of the double Higgs production cross sections. In particular, the H —>■ hh 
resonant enhancement was the main focus of our studies, involving the lightest, h, and the 
heaviest, H, of the neutral Higgs bosons of the MSSM, in the kinematic regime m# > 2rrih- 
This dynamics can for example occur in the following reactions: gg — > hh in the hadronic 
case and e + e~ — > hhZ in the leptonic one, but only at low tan/3. These two processes 
proceed via intermediate stages of the form gg —> H and e + e~ — > HZ, respectively, followed 
by the decay H — > hh. Thus, they in principle allow one to determine the strength of the 
Hhh vertex involved, \nhh-, in turn constraining the form of the MSSM Higgs potential itself. 
The signature considered was hh — > bbbb, as the h — > bb decay rate is always dominant. 

We have found that several kinematic cuts can be exploited in order to enhance the 
signal-to-background rate to level of high significance, particularly at the e + e~ machine. At 
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the pp accelerator, in fact, the selection of the signal is made much harder by the presence 
of an enormous background in 4b final states due to pure QCD. In parton level studies, 
based on the exact calculation of LO scattering amplitudes of both signals and backgrounds 
(without any showering and hadronisation effects but with detector acceptances), we have 
found very encouraging results. At a LC, the double Higgs signal can be studied in an 
essentially background free environment. At the LHC, the signal and the QCD background 
are in the end at the same level with detectable but not very large cross sections. 

Earlier full simulations performed for the e + e~ case had already indicated that a more so- 
phisticated treatment of both signal and backgrounds, including fragmentation/hadronisation 
and full detector effects, should not spoil the results seen at the parton level. For the LHC, 
our preliminary studies of gg — > H — > hh — > bbbb in presence of the gg — > hh — > bbbb 
continuum (and relative interferences) also point to the feasibility of the signal selection, 
after realistic detector simulation and event reconstruction. As for double h production in 
the continuum, although not very useful for Higgs self-coupling measurements, this seems 
a promising channel, if not to discover the lightest MSSM Higgs boson certainly to study 
its properties and those of the Higgs sector in general (because of the large production and 
decay rates at high tan f3 and its sensitivity to such a parameter), as shown from novel simu- 
lations also presented in this study. (The discovery potential of this mode will eventually be 
addressed in Ref. 0.) Despite lacking a full background analysis in the LHC case, we have 
no reason to believe that a comparable degree of suppression of background events seen at 
parton level cannot be achieved also at hadron level. Progress in this respect is currently 
being made ||. 
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